Solid state quantum memory using the 31 P nuclear spin 
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The transfer of information between different physical 
forms is a central theme in communication and computa- 
tion, for example between processing entities and mem- 
ory. Nowhere is this more crucial than in quantum com- 
putation PQ, where great effort must be taken to pro- 
tect the integrity of a fragile quantum bit (qubit) [2]. 
However, transfer of quantum information is particularly 
challenging, as the process must remain coherent at all 
times to preserve the quantum nature of the informa- 
tion [3J. Here we demonstrate the coherent transfer of a 
superposition state in an electron spin 'processing' qubit 
to a nuclear spin 'memory' qubit, using a combination 
of microwave and radiofrequency pulses applied to 31 P 
donors in an isotopically pure 28 Si crystal [4j [5]. The 
state is left in the nuclear spin on a timescale that is 
long compared with the electron decoherence time and 
then coherently transferred back to the electron spin, 
thus demonstrating the 31 P nuclear spin as a solid-state 
quantum memory. The overall store/readout fidelity is 
about 90%, attributed to imperfect rotations which can 
be improved through the use of composite pulses [6 . The 
coherence lifetime of the quantum memory element at 
5.5 K exceeds one second. 

Classically, transfer of information can include a copy- 
ing step, facilitating the identification and correction of 
errors. However, the no-cloning theorem limits the abil- 
ity to faithfully copy quantum states across different de- 
grees of freedom [7]; thus error correction becomes more 
challenging than for classical information and the trans- 
fer of information must take place directly. Experimental 
demonstrations of such transfer include moving a trapped 
ion qubit in and out of a decoherence-free subspace for 
storage purposes [8] and optical measurements of NV cen- 
tres in diamond [9 . 

Nuclear spins are known to benefit from long coher- 
ence times compared to electron spins, but are slow to 
manipulate and suffer from weak thermal polarisation. 
A powerful model for quantum computation is thus one 
in which electron spins are used for processing and read- 
out while nuclear spins are used for storage. The storage 
element can be a single, well-defined nuclear spin, or per- 
haps a bath of nearby nuclear spins [10 . 31 P donors in sil- 
icon provide an ideal combination of long-lived spin- 1/2 
electron [11] and nuclear spins [12] , with the additional 



advantage of integration with existing technologies [4 
and the possibility of single spin detection by electrical 
measurement [I3j [HJ US]- Direct measurement of the 
31 P nuclear spin by NMR has only been possible at very 
high doping levels (e.g. near the metal insulator tran- 
sition [16]). Instead, electron- nuclear double resonance 
(ENDOR) can be used to excite both the electron and 
nuclear spin associated with the donor site, and measure 
the nuclear spin via the electron [17]. This was recently 
used to measure the nuclear spin-lattice relaxation time 
Ti n , which was found to follow the electron relaxation 
time Ti e over the range 6 to 12 K with the relationship 
Ti n ^ 250Ti e [5] [12]. The suitability of the nuclear spin 
as a quantum memory element depends more critically 
on the nuclear coherence time T2 n , the measurement of 
which has now been made possible through the storage 
procedure described here: by varying the storage time 
and observing the amplitude of the recovered electron 
coherence. 

Figure [TjB) shows the coherence transfer scheme used 
for the write process from a processing qubit represented 
by an electron spin degree of freedom, to a memory qubit 
residing in a nuclear spin degree of freedom. Each tt pulse 
is equivalent to a controlled-NOT gate [18] (with some 
additional phase which can be ignored) such that the 
pair of 7r pulses constitute a SWAP gate. The scheme as- 
sumes that all pulses are on-resonance and have sufficient 
bandwidth to completely excite an individual transition. 
A read operation is performed by applying the reverse 
sequence to bring the coherent state back to the electron 
spin qubit. Although the phase relationship betweeen 
the microwave and rf pulses must be constant through- 
out this process, any phase difference is cancelled out 
over the course of the write-read process. In practice, 
this means the microwave and rf sources need not be 
phase-locked, but must have high phase stability. This is 
illustrated in calculations following the evolution of the 
density matrix, provided in the Supplementary Material. 

Although the electron spin qubit can be prepared in 
a state of high purity using experimentally accessible 
magnetic fields and temperatures, the small nuclear Zee- 
man energy results in the nuclear spin being initially in 
a highly mixed thermal state. However, for the purposes 
of this quantum memory scheme it is not necessary to 
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FIG. 1: The level structure of the coupled electron and 
nuclear spins and scheme for the transfer of a logical 
qubit within the two physical spin qubits. (A) The four 
level system may be manipulated by resonant microwave and 
radiofrequency (rf) radiation. In our experiments the logical 
electron spin 'processing' qubit is represented by states |1) and 
1 2), whose state can be transferred to a nuclear spin 'memory' 
qubit represented by states |2) and |4). State |3) is never 
addressed at any point and can be ignored. (B) An electron 
spin coherence between states |1) and |2) is transferred to the 
nuclear spin qubit by an rf tt pulse followed by a microwave 
7r pulse. Both pulses must fully excite the transition, and 
be short compared with the electron and nuclear coherence 
times. The reverse process is used to transfer the nuclear 
coherence back to the electron. 
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FIG. 2: Coherent storage of an electron spin state in a 
nuclear spin state, using 31 P-doped 28 Si-enriched sili- 
con single crystal. A) An electron spin coherence is stored 
in the nuclear spin for 2r n ~ 50 ms, at 7.2 K. The recovered 
electron spin echo is of comparable intensity to that obtained 
at the beginning of the sequence, even though the electron 
spin coherence time T2 e here is about 5 ms. The lifetime of 
the stored state is limited instead by the nuclear decoherence 
time T2n, which can be measured directly by varying r n . B) 
The recovered echo intensity was measured a function of the 
storage time at 5.5 K while applying a dynamic decoupling se- 
quence (CPMG) to the nuclear spin, yielding a exceeding 
1 second. 



perform any pre-cooling of the nuclear spin resource [34] . 

The above model is sufficient given a single electron- 
nuclear spin pair, or a homogenous ensemble. However, 
in the experiment described here, we must consider the 
effects of inhomogeneous broadening across the ensem- 
ble of spins being manipulated. The effect of inhomo- 
geneous broadening is to leave some electron (nuclear) 
spins detuned from the applied microwave (rf) radiation, 
by S e (S n ). In a suitable rotating reference frame, electron 
(nuclear) spin coherence will thus acquire an additional 
phase at a rate S e (5 n ), while double quantum coher- 
ences will acquire phase at a rate S e + S n . Thus, inho- 
mogeneous broadening requires the application of care- 
fully placed refocusing pulses to bring all spin packets 
into focus at key points during the transfer process. In 
the experiment described here, 7r/5 e ~ 2 /xs and 7r/5 n ~ 
100 /is. 

Figure [2] shows the practical implementation of a pro- 
tocol that generates a coherent electron spin state, stores 
it in a state of the nuclear spin for some time, and then 
recovers it to the electron state for readout again. The 
coherence is first generated by a microwave tt/2 pulse of 
a chosen phase representing our bit of quantum in- 
formation. A free induction decay (FID), the reversible 
dephasing of the ensemble, follows this pulse. We apply 
a refocusing microwave tt pulse at time r e to initiate a re- 
vival in the electron spin coherence. The subsequent rf tt 
pulse transfers the coherence from the electron to a dou- 



ble quantum coherence of entangled electron-nuclear spin 
states. During this period the phase 5 e T e , acquired be- 
fore the microwave refocusing pulse, continues to unwind 
so that when the final step of the transfer, a microwave 
tt pulse, is applied the effect of the inhomogeneous elec- 
tron spin packets has been completely refocused. The 
quantum information that was generated by the first mi- 
crowave tt/2 pulse now resides entirely in the state of the 
nucleus. 

This information may be stored in the nuclear state for 
some extended period so the effects of inhomogeneities on 
the phase of the nuclear state become appreciable and a 
preparatory rf refocusing pulse must be applied before 
the information can be recovered. During the nuclear 
spin echo, the coherence is transferred back to the elec- 
tron state with a microwave tt pulse followed by an rf 
tt pulse. We apply one further microwave tt pulse to 
stimulate an electron spin echo representing the readout 
event. The lower right panel of Figure [2] shows the real 
(red) and imaginary (black) parts of this echo for dif- 
ferent initial phases cp, demonstrating that the encoded 
phase is recovered through the storage-recovery process, 
as required for an effective quantum memory element. 

The storage time is limited only by the nuclear deco- 
herence time T2 n , which is in turn limited to 2Ti e when 
there is a significant hyperfine interaction (A ^> 1/Ti e ) 
between the electron and nuclear spin and in the low- 
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FIG. 3: Observing the nuclear spin coherence during 
the storage process. The phase of the initial electron su- 
perposition state is determined by the phase of initial tt/2 mi- 
crowave excitation pulse, which we can control. This state is 
then transferred to the nuclear spin using the scheme outlined 
in Figure [I] The nuclear spin coherence is read using a pro- 
cess similar to a Ramsey fringe measurement: an rf tt/2 pulse 
converts nuclear coherence to nuclear polarisation, which is 
then detected via an electron spin echo measurement selec- 
tive to one nuclear spin state. The correlation of the phase of 
the nuclear spin echo to the phase of the original electron spin 
superposition confirms the coherent nature of the transfer for 
electron to nuclear spin. 



field/high-temperature limit (see Supplementary Mate- 
rial); Ti e becomes very long (e.g. hours) at low temper- 
atures [17]. A direct measurement of T 2n in anything 
other than highly-doped Si:P has been impossible by 
traditional NMR means, but our write/read procedure 
provides a method for performing this measurement by 
increasing the storage time T st0 re and observing the re- 
sulting decay in the recovered electron coherence. T 2n ob- 
tained in this way indeed follows 2Ti e approximately over 
the range 9 to 12 K as expected, though at lower tem- 
peratures an additional T 2n process appears to play a 
role, yielding a limit of about 65 ms. A leading can- 
didate for this additional process is slowly fluctuating 
fields, the effect of which may be mitigated by dynami- 
cally decoupling the system [19] [20] . By applying a Carr- 
Purcell-Meiboom-Gill (CPMG) decoupling sequence [21] 
at a 1 kHz repetition rate to the nuclear spin during the 
storage period, we were able to obtain much longer deco- 
herence times than for a simple Hahn echo measurement, 
rising to 1.75 seconds at 5.5 K, as shown in Figure [2^. 

Under optimised conditions, T 2e is limited only by 
magnetic dipole-dipole interactions, and values between 4 
and 6.5 ms have been measured in the samples used here, 
varying according to the donor spin concentration [TT] . 
Using the nuclear degree of freedom, we have achieved 
storage times several orders of magnitude longer than 
T 2e . 



The removal, or substantial detuning of any of the rf 
pulses in the sequence destroys the recovered echo, con- 
firming the importance of the transfer to the nuclear spin 
and providing evidence that the stored quantum informa- 
tion does indeed reside in the nuclear state. To go further 
we require a tool permitting introspection of the state of 
the nuclear spin during the storage period. We therefore 
applied a sequence to (destructively) probe the nuclear 
coherence via the electron state, as shown in the upper 
panel of Figure [3] The early part of the sequence is as 
described above: an electron spin coherence is stored in 
the state of the nucleus. When we would like to observe 
the state of the nucleus, we apply an rf tt/2 pulse to 
convert the nuclear coherence into a nuclear polarisation 
(in the spirit of a Ramsey fringe experiment). A short 
electron spin echo sequence, selective in one nuclear sub- 
space, then reveals the population of the nuclear level. 

This sequence can be performed at any time; the lower 
panel of Figure [3] shows the result of observing the state 
of the nucleus at a range of times for different starting 
phases (/?, revealing the nuclear spin echo following the 
rf refocusing pulse. The centre of the rf frequency was 
intentionally moved off-resonance to produce oscillations 
in the nuclear echo to aid the identification of the phase 
of the nuclear coherence. The fact that the phase of the 
nuclear spin echo follows the phase of the original mi- 
crowave tt/2 pulse confirms that the information transfer 
process has remained coherent [35] . 

To demonstrate the generality of the storage sequence 
described here, we applied it to a wider set of initial 
states, in particular the ±X, ±Y, ±Z and Identity basis 
states and performed density matrix tomography by com- 
paring the original states with those recovered after the 
write-read process (see Supporting Information for full 
details). The results are summarised in Figure [4j and 
show fidelities of approximately 0.90, where the fidelity 
between initial (pseudo) pure state po and recovered state 
pi is denned as F = pi |^), where po = We 
attribute the reduced fidelity to a ^ 5% error in each 
of the seven microwave and rf pulses applied over the 
course of the sequence, which is entirely consistent with 
previous measurements of pulse fidelities [22]. Such er- 
rors are mostly systematic, and may be corrected through 
the application of composite pulses, as previously demon- 
strated in both EPR and NMR [6j[23]. By replacing some 
of the microwave pulses with BB1 composite pulses we 
were able to improve the overall fidelity to approximately 
0.97 and further improvements are to be expected with 
greater control of the rf pulse phases. 

As the experimental challenges of quantum informa- 
tion processing have become better understood, the im- 
portance of hybrid quantum systems in models for quan- 
tum information has emerged [24j [25j [26j [27]. The 
approach described here demonstrates the advantages 
of such hierarchical models and has a broad applica- 
bility in systems where there is a substantial asymme- 



4 



+X +Y +Z Identity 




Fidelity: 0.905 0.887 0.881 1.000 



FIG. 4: Density matrix tomography for original and 
recovered states. Pseudopure states +X, +Y and +Z, and 
the Identity were prepared in the electron spin qubit and mea- 
sured (first row) . These states were then stored in the nuclear 
spin degree of freedom and then returned to the electron spin 
and measured (second row). Tomography was performed by 
measuring the qubit in the a x , a y and a z bases. The fidelity of 
the quantum memory was obtained by comparing the initial 
and recovered density matrices. 



try in relaxation times. Storage can be driven glob- 
ally, as shown here, or locally, using EPR gates [28] or 
Stark tuning [4 j. Furthermore, our protocol for faith- 
fully transferring a coherent electron spin state to the nu- 
clear spin offers a route to projective measurements of the 
qubit state through proposed spectrally sensitive single- 
spin-detection methodologies such as STM-detected or 
electrically-detected EPR [29]. 
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METHODS 



Si:P consists of an electron spin S=l/2 (g = 1.9987) coupled to the nuclear spin / = 1/2 of 31 P through a hyperfine 
coupling A = 117 MHz [17], and is described by an isotropic spin Hamiltonian (in angular frequency units): 

H =u; e S z -u; I I z + A.S.T, (1) 

where oj e = g/3Bo/h and uoi = giP n Bo/h are the electron and nuclear Zeeman frequencies, g and gi are the electron 
and nuclear g-factors, (3 and (3 n are the Bohr and nuclear magnetons, H is Planck's constant and Bo is the magnetic 
field applied along z-axis in the laboratory frame. The X-band EPR signal comprises two lines (one for each nuclear 
spin projection Mi = ±1/2). Our experiments were performed on the high-field line of the EPR doublet corresponding 
to Mi = -1/2. 

Single crystal samples were used, as epilayers of 28 Si have a biaxial residual stress that broadens the 31 P ENDOR 
line and makes it difficult to fully excite. 28 Si-enriched single crystals with a residual 29 Si concentration of 800 ppm 
were produced by decomposing isotopically enriched silane in a recirculating reactor to produce poly- Si rods, followed 
by floating zone crystallization [30]. To reduce spin-spin coupling effects, the phosphorus concentration was reduced 
from an initial value of near 1 • 10 15 cm -3 to 2-5- 10 14 cm -3 by five passes of zone refining followed by floating zone 
crystallisation. 

Pulsed EPR experiments were performed using an X-band (9-10 GHz) Bruker EPR spectrometer (Elexsys 580) 
equipped with a low temperature helium-flow cryostat (Oxford CF935). The temperature was controlled with a 
precision greater than 0.05 K using calibrated temperature sensors (Lakeshore Cernox CX-1050-SD) and an Oxford 
ITC503 temperature controller. 

For most measurements, microwave pulses for tt/2 and tt rotations of the electron spin were set to 700 and 1400 ns, 
and no travelling wave tube (TWT) amplifier was used. For CPMG and BB1 experiments, an Amplifier Research 
20W solid state CW amplifier was used, with tt/2 and tt pulses 80 and 160 ns respectively. RF pulses of 20 ps were 
used for tt rotations of the 31 P nuclear spins. During CPMG, up to 1000 refocusing pulses were applied during a 
single sequence. 



SUPPLEMENTARY INFORMATION 



A. Effect of radiofrequency and microwave phases 



The chosen basis is: 



1 1 
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where S represents the electron donor spin, and I the 31 P nuclear spin. All pulses are assumed to be selective on a 
particular electron or nuclear spin transition, as illustrated in Figure [I] The phase of the initial tt/2 microwave pulse, 
y? e , determines the phase of the initial electron spin coherence. All other microwave pulses have phase (p mw , while 
that of the rf pulses is <p r f. 

The initial spin density matrix, neglecting any nuclear spin polarisation, is proportional to (I + (3S Z ), where I is 
the Identity matrix and (3 = — gfl ^° . This can be rewritten in pseudopure state form, neglecting the majority of the 
Identity component and omitting the constant factor (3: 



Po = (S z +I/2)/2 = Pth 



( 1/2 \ 



1/2 

V o o o o J 



(3) 



After the initial (coherence-generating) tt/2 microwave pulse: 



Pi 



( 1/4 


exp(-i<p e )/4 





o\ 


exp(z<p e )/4 


1/4 














1/2 





V o 











(4) 



7 



The next two pulses, ttrf followed by ir mw , transfer this coherence to the nuclear spin: 
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The coherences here decay with characteristic time T2 n , which is typically much longer than T^ Q . Upon applying 
the reverse of the transfer sequence above (7r mw followed by 7r r j), the electron coherence is revived: 
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(6) 



Thus the relative phase of the microwave and rf sources is cancelled out, though both must remain stable over the 
course of the experiment. 



B. Electron and nuclear spin relaxation 



Electron relaxation of the system can be modeled using a standard master equation in Lindblad form. In order to 
represent processes that take the system to thermal equilibrium, both raising and lowering terms are included. 

p = -^(pS~S+ + S'S+p - 2S+S-) - ^^(pS+S- + S+S'p - 2S~S + ) - i[H, p] (7) 

where 7 is the relaxation rate, S + and S~ are the electron spin raising and lowering operators (S ± = S x db iS y ), and 
(3 relates the electron Zeeman splitting to /c#T, as defined above. This can be simplified by transforming into the 
rotating frame of the Hamiltonian, taking an Ising approximation (Ho = oo e S z — ooil z + A-S z -I z ). In this frame, H 
goes to 0, and we are left only with the relaxation part of Eq. (7), with 5 + and S~ transformed into the rotating 
frame. Neglecting direct nuclear relaxation (Ti n — > 00) and in the high temperature limit this yields: 
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(8) 



The electron relaxation rate (1/Ti e ) can be ascertained by observing the appropriate density matrix elements: 
Al,i + P3,3 

Taking p e = p lyl + p 3j3 then, p e = - j(p e - 1/2). Solving this gives: 



T T 1 1 

- ^ (Pl,l + PS,S) + 2^2 + P4,4) = - ^ (Pl,l + PsA + 2 ( X ~ ^ ~ ^ 3 ' 3 ) 



Pe = Pe, e" 7 ' + 1/2 

Hence, electron relaxation follows e~ 7t and the electron relaxation time, T\ e 



(9) 



The nuclear coherence is given by p nn = p 3 ^i + P4,2- Extracting these terms from Eq.|8]yields two coupled differential 
equations: 



P3,l 
P4,2 



1 -e 

— iAt 2 



P3,l 
P4,2 



(10) 



It is straightforward to get rid of the time dependence in the 2x2 matrix in this equation by making a time 
dependent unitary transformation and solving for new variables p' 3 x and p' A 2 - 



P'3,1 

P4,2 



u 



Ps,i 

P4,2 



-iAt /I 





JAt/2 



PS,1 
P4,2 



(ii) 
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FIG. 5: Pulse sequences applied to prepare and measure electron spin states. 



Following this transformation solving the pair of differential equations is a simple eigenvalue problem. In the exper- 
iments A = 117 MHz and 7 ranges from 1 kHz to less than 1 Hz (as a function of temperature), hence we can take 
the limit A ^> 7. In this case, both characteristic eigenvalues have a real part of —7/2, and therefore any nuclear 



coherence decays with this rate. Thus T2 n 



2Ti e as experimentally observed. 



C. Density matrix tomography 



In this section we describe 1) the preparation and the tomography of the pseudopure initial electron spin states 
±X, ±y, and ±Z, and the Identity, 2) the tomography of the state after transfer to and from the nuclear spin and 
3) the measure of fidelity between the starting and recovered states. 

We define state ±X as (±a x + I)/2, and similarly for Y and Z. Figure [5] shows the full set of pulse sequences 
required for the state preparation and detection. The starting (thermal) state is +Z, which thus requires no preparation 
pulse. — Z is obtained by applying an inversion tt pulse, while ±X and ±Y are obtained through a n/2 pulse of the 
appropriate phase. The state I is obtained by applying an inversion tt pulse and waiting some time, T = (ln2)Ti e . 
This is long enough to ensure complete decoherence of the electron spin (off diagonal elements go to zero), while 
corresponding to the precise point during the relaxation process at which the electron spin populations are equal. 

Measurement is performed in the a x and a y bases by generating an electron spin echo and observing both in-phase 
and quadrature components. A measurement in the a z basis can be performed some short time (t < T ie ) later in 
the pulse sequence by applying a tt/2 pulse, followed by a tt pulse, and then observing the resulting echo. This 
measurement operation is applied to both the starting states and those recovered after the end of the write-read 
process. 

Echo traces for the seven states are shown in Figure [6] for both the starting and recovered states. Each corresponds 
to a measurement in the cr x ,y,z bases. 

The integrated areas of the electron spin echoes, A Xj y iZ , are used to extract the components of <J x ,y,z m the spin 
density matrix. We assume the starting electron spin state is (pseudo)pure and can thus normalise the areas to extract 
a density matrix of the starting electron spin state: 



A X (J X + Aydy + A Z (J Z . 

2y A* + A^ + A 2 Z 



(12) 
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+X -X +Y -Y +Z -Z Identity 




FIG. 6: Measurements in the cr x ,y,z bases of the electron spin state before and after storage in the nuclear 
spin. Electron spin echoes (a.u.) are obtained using the pulse sequences shown in Figure [5] The echoes for a x and a y occur 
simultaneously, while that for a z (which occurs at a later time) is superimposed here for clarity. 



+X +Y +Z Identity -X -Y -Z 




Fidelity: 0.905 0.887 0.881 1.000 0.905 0.895 0.885 



FIG. 7: Density matrices for the initial and recovered states. 

We make no such assumptions about the purity of the recovered electron spin state, and normalise the integrated 
areas of the recovered spin echoes using the areas obtained from the starting state. Density matrices for the initial 
and recovered state are shown in Figure [7] 

One common measure of the difference between two quantum states is the fidelity [31] : 

F(po,pi) = Tr(y/y/pipoVpi) (13) 

Here, we use a more aggressive measure of fidelity, F' = F 2 , corresponding to the overlap of a pure state and an 
arbitrary density matrix (rather than its square root). Thus, if our initial density matrix po = the fidelity 

measure we use is: 



F' = {^\ Pl \i>). 



(14) 



